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ABSTRACT 

We study the relationship between the metallicity of gamma-ray burst (GRB) progenitors and the probability 
distribution function (PDF) of GRB host galaxies as a function of luminosity using cosmological hydrodynamic 
simulations of galaxy formation. We impose a maximum limit to the gas metallicity in which GRBs can occur, 
and examine how the predicted luminosity PDF of GRB host galaxies changes in the simulation. We perform 
the Kolmogorov-Smirnov test, and show that the result from our simulation agrees with the observed luminosity 
PDF of core-collapse supernovae (SNe) host galaxies when we assume that the core-collapse SNe trace star 
formation. When we assume that GRBs occur only in a low-metallicity environment with Z < O.IZq, GRBs 
occur in lower luminosity galaxies, and the simulated luminosity PDF becomes quantitatively consistent with 
the observed one. The observational bias against the host galaxies of optically dark GRBs owing to dust 
extinction may be another reason for the lower luminosities of GRB host galaxies, but the observed luminosity 
PDF of GRB host galaxies cannot be reproduced solely by the dust bias in our simulation. 
Subject headings: gamma rays: bursts — galaxies: formation — method: numerical 



1. INTRODUCTION 

Long gamma-ray bursts (GRBs) are the brightest astronom- 
ical transient events, and they are important laboratories of 
high energy astrophysics in extreme conditions, as well as the 
tools to probe the high-redshift universe. For example, peo- 
ple have used GRBs to estimate the cosmi c star formation rat e 
(SFR) density at very high redshifts (e.g.. lKistler et al.ll2009l) . 
The observational association of some of t he long GRBs with 
energetic Type Ic supernovae (SNe, e.g.. iHjorth et al.l [20031: 
IStanek et al.l 120031) is considered to be the evidence that at 
least some of the GRBs originate from core-collapses of very 
massive stars. Zhang et al. (2009) applied multiple criteria 
to identify GRBs whose progenitors are related to core col- 
lapses of massive stars (the so-called "Type 11" GRBs), and 
found that they mostly correspond to the traditional 'long' 
GRB population. Hereafter 'GRB' means Type II GRBs, and 
'SN' means core-collapse SN, unless otherwise stated. 

Although GRBs may originate from the core-collapses of 
very massive stars, the occurrence rate of GRBs is much lower 
than that of normal SNe, and the conditions required for a 
GRB to occur from a SN still remain as one of the most out- 
standing questions in current astrophysics. A better under- 
standing of required conditions for a GRB to occur would 
constrain the physical mechanisms of GRBs, and allow us 
to predict the GRB detection rate for the future observations 
more reliably. 



Some theoretical studies on the origin of GRBs using stel- 
lar evolution models suggest that a low metallicity envi- 
ronment may be a necessary condition for a GRB to oc- 
cur (e.g., iMacFadyen & Wooslevll99llYoon & Langed2005l: 
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IWooslev" &Hegeiil2006ir It has also been suggested from 
the observations that the metallicity distribution of GRB host 
galaxies at redshift z < 0.25 is significantly biased towards 
low metallicities compared to the expectati on when GRBs 
are unbiased tracer s of star formation (Stanek et al.l 120061 
iModjaz et al.l l2008h . However, reliable spectroscopic esti- 
mates of metallic ities are available on ly for galaxies at low 
redshifts (z < 0.5: ISavaglio et al]l2009h while the majority of 
GRBs occur at higher redshift. 

Furthermore, some observations suggest that the GRB host 
galaxies are system atically fainter th a n those of the core- 
collapse SNe (Le Floc'h et all llOOl iFruchter et all l2006l 
hereafter F06), indicating that the GRBs may preferentially 
occur in low metallicity environment, because fainter and 
lower mass galaxies generally have lower metallicities. These 
interpretations have also been supported by other theoretical 
studies usin g the mo dels of galaxy formatio n and evolution 
(iNuza etal.li 2007: La pi et al.ll2008t ICampisi et al. 2009). 

It is also reported that the GRB host galaxies at z > 2 
have larger Lya equivalent widths comp ared to general star- 
forming galaxies at similar redshifts ( Jakobsson et al.l 120051: 
|Fynbo et al. 2002, 2003), which is possibly a result of the 
stronger ionizing flux emitted from low-meta llicity stellar 
population in GRB host galaxies (Niin o et ani2009h . Al- 
though the Lya emission property of GRB host galaxies could 
be used as a metallicity indicator in the studies of GRB pro- 
genitors, our current understanding of Lya transfer in the in- 
terstellar medium (ISM) is still inadequate to draw a robust 
conclusion ( Niino et a l. 2009), and the samples of GRB host 
galaxies with Lya detections is currently very small. 

The differences between the observed probability distribu- 
tion functions (PDF) of GRB and SN host galaxies at z ~ 1 
as a function of luminosity (hereafter 'luminosity PDF') or 
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size is consistent with the hypothesis that GRBs occur prefer- 
entially in low-metallicity environment (F06). However, it is 
not straightforward to connect the difference in the luminos- 
ity and/or size of host galaxies to the metallicity difference 
of GRB and SN progenitors. In fact, some studies of GRB 
host properties claim that GRBs can be produced in higher 
metallicity envi ronments than suggested by the stellar evo- 
lution models ( Wolf & Podsiadlowskil 120071 : iKocevski et alj 
120091). Ho wever, their results suffer from some uncertainties 
(see § 14. 2b . and more tests are required to establish the con- 
nection between the metallicity dependence of GRBs and the 
luminosity PDF of GRB host galaxies. 

In this paper, we investigate the relationship between the 
low metallicity preference of GRBs and the luminosity of 
GRB host galaxies at z ^ 1, where relatively large sample of 
observed GRB/SN host galaxies is available, using cosmolog- 
ical hydrodynamic simulations. We compare the predictions 
of our simulations with the observations to test if the low 
metallicity preference of GRBs predicted by the stellar evo- 
lution models is quantitatively consistent with the observed 
luminosity PDF of GRB host galaxies. Several studies on 
this topic using galaxy formation models have already been 
carried out (Nuza et al. 2007; Lapi et al. 2008; Campisi et aD 
120091) . However, quantitative comparisons between the pre- 
dicted luminosity PDF of GRB host galaxies and the observed 
luminosity PDF have not been performed yet. 

This paper is organized as follows. In §|2l we briefly de- 
scribe our simulation code and show some relevant properties 
of simulated galaxies, such as the luminosity function (LF), 
dust extinction, and luminosity-metallicity relationship. In 
§[3] we describe the modeling of GRB/SN event rate in the 
simulated galaxies, and compare the resulting luminosity PDF 
of the GRB/SN host galaxies to the observations. In §21 we 
discuss the effect of possible observational bias by the opti- 
cally dark GRBs on the luminosity PDF of GRB host galax- 
ies. Then we compare our results with previous studies. Our 
conclusions are summarized in §|5] 

2. NUMERICAL METHOD AND BASIC DATA 
2.1. Simulations and Galaxy Identification 

We use the modified version of the tree-particle-mesh 
smoothed particle hydro dynamics (SPH ) code GADGET-3 
(originally described in SpringelM 2005h . In this code, the 
SPH calculation is performed base d on the en tropy conserva- 
tive formulation (ISpringel & Hernquistll2002h . Our conven- 
tional code includes radiative cooUng by H, He, and metals 
jChoi & Na gamine 2009) , heating by a unif orm UV back- 
ground of a modified Haardt & Madaul ( 11996 ) spectrum, star 
formation, supernova feedback, phenomenological model for 
galactic winds, and a sub-resolution model of multiphase ISM 
and star formation (Springel & Hernquist 20031). 

In this paper, we adopt the following fiducial cosmology 
which is consistent with the latest WMAP result: ri,„ = 0.26, 
VL^ = 0.74, = 0.044, h = 0.72, n, = 0.96, and erg = 0.80 
dKomatsu et al.ll2009l 12010 1). We use two runs with different 
box sizes and resolution: N400L34 and N400L100 run. See 
Table[T]for the simulation parameters. 

We identify galaxies in the simulation at z = 1.0 using a 
simpUfied variant of the S UBFIND algorithm (ISpringel et al.l 
llOOltlChoi & Nagamineir2009) . In more detail, the code first 
computes a smoothed baryonic density field to identify can- 
didate galaxies with high density peaks. The full extent of 
these galaxies are found by adding gas and star particles to 



TABLE 1 
Simulation Parameters 



Run 


Box-size 




"iDM 


'«ga; 




€ 


N4()()L34 


33.75 


2 X 400^ 


3.60 X 10' 


7.33 X 


10* 


3.375 


N400L100 


100.0 


2 X 400^ 


9.12 X 10** 


1.91 X 


10« 


6.45 



Note. — Simulations employed in the present work. The box-size is 
given in units of /r' Mpc, Np is the pailicle number of dark matter and gas 
(hence X 2), moM and mgas are the masses of dark matter and gas parti- 
cles in units of /!~'Mq, respectively, and e is the comoving gravitational 
softening length in units of A"' kpc. The value of e is a measure of spatial 
resolution. 

the galaxies in the order of declining density. If all Nmin near- 
est neighbor particles have lower densities, this group of par- 
ticles is considered as a new galaxy. Here, N„,i„ is the min- 
imum number of gas and star particles that constitute one 
isolated galaxies. In this paper we set N,„i„ = 32. If there 
is a denser neighbor, the particle is attached to the galaxy 
to which its nearest denser neighbor already belongs to. If 
two nearest neighbors belong to different galaxies and one of 
them has less than N,„i„ particles, then the two galaxies are 
merged. If two nearest neighbors belong to different galax- 
ies and both of them has more than A',,,,,, particles, the parti- 
cles are attached to the larger galaxy, leaving the other galaxy 
intact. In addition, the gas particles in galaxies should be 
denser than O.Olp,/,, where p,i, is the star formation density 
threshold (S pringel & He rnquist 2003). In the current sim- 
ulations, we use the SF threshold density of n,/, = 0.6 cm"^ 
dChoi & Nagaminel2010l) . 

2.2. Distribution of Gas and Stars 

We first examine the distribution of gas and stars in some 
simulated galaxies. Figure [T] is an example of a typical faint 
GRB host galaxy with Muv = -19.0, while Figure|2]is a more 
luminous galaxy with Muv = -20.3. In the lower panels of 
Figures [T] and |2l the distribution of gas with high- and low- 
metallicity is shown separately. In both examples, the high- 
metallicity gas is concentrated at the center of its host galaxy, 
and the low-metallicity gas is distributed more broadly. 

It is reported that the observed GRBs are primarily located 
in the brightest regions of galaxies compared to SNe (F06), 
and the concentration of high-metallicity gas at the center of 
galaxies in our simulation may seem to contradict with the 
low-metallicity preference of GRBs. However, if the GRBs 
occ ur in the young star c lusters that are not heavily enriched 
yet (Larsso n et al.ll20()7l) . then their occurrence in bright re- 
gions would be unrelated to metallicity. We also see from this 
figure that our simulations do not have adequate resolution to 
resolve the disk structure with spiral arms, therefore it is dif- 
ficult to discuss the spatial distribution of GRB sites within a 
single galaxy using our current cosmological simulations. 

2.3. Dust Extinction Model 

We compute the spectra and luminosities of simulated 
galaxi es by applying the GA LAXEV population synthesis 
code (Bruzual & Charlol 12003) to each constituent star par- 
ticles and employing a simple model of dust extinction de- 
scribed below. To estimate the degree of extinction in each 
simulated galaxy, we assume that Eg-v is proportional to the 
metal mass column density of the ISM: (Jism x Zism, where 
fJiSM and ZisM are the gas column density and metallicity 
of each galaxy, respectively. The proportionaHty constant 
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Fig. 1. — Distribution of gas {top left) and stars (top right) in an example 
galaxy from the N400L34 run. This galaxy has Muv = -19.0, Mg = -19.3, 
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M^, = 3.5 X 10'Mq, and 

2.0Zgaiaxy and Z < O.lZgaiaxy are shown in the bottom left panel and the bot- 
tom right panel, respectively. 
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Fig. 2. — Same as Fig.[T] but for a more luminous galaxy. This galaxy 
is from the N400L100 simulation, and has the following properties: Muv = 
-20.3, Mg =-21.8, =4.5 X 10"Mo, andZgai^xy = 1.2Zq. 



is fixed so that it agrees with the empirical estimate of gas- 
to-dust ratio in th e Milky Way {Nu/Eb-v = 5.8 x 10^' cm"^; 
iBohlin et al.ll 19781) when Zjsm = Z©: 



/h ctism Zism 



5.8 X 1021 



(1) 



where /h = 0.75 is the hydrogen mass fraction and nip 
is the proton mass. This formulation also agrees with 
the gas-to-dust ratio in the Small Magellanic Cloud (SMC) 
Nh/Eb-v ~ 10^^ cm-2 (e.g.. iBouchet et al.lll985l: .FitzpatricS 
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Fig. 3. — Extinction vs. metallicity for the simulated galaxies. The solid 
(red) and dashed (blue) contours are for the N400L34 and N400L100 runs, 
respectively. The contours represent ng^\/ dV dZdEg-v = 0.01, 0.1,10, 100, 
and 1000 [/i^Mpc^^n^ag"']- 

IT981 iTumlinson et al.l l2002h . given that the metallicity in 
SMC is - 1O"'Z0. We use the ICalzettil (11997b extinction 
curve to calculate the extinction at each wavelength. 

To estimate ctism for each simulated galaxy, we assume 
that the ISM in each simulated galaxy follows the Hernquis^ 
([1990) mass profile with a total gas mass Mism of the galaxy. 
Then the mean ctism within a radius r is 



CTlSM = 



MHeraquist(< r) _ MiSM 



87r(r+r,.)2' 



(2) 



where the characteristic radius of the Hernquist mass profile 
is fixed such that MHemquist(< ^smooth,o) = '«part,o- Here, / 

.smooth.O 

and mpait.o are the SPH smoothing length and the mass of the 
central gas particle with the highest gas density in the galaxy. 
For simplicity, we further assume ctism = A#iSM/87rr^, which 
gives the value in the central region of the galaxy (r ^ rc), 
and set Zism to the metallicity of^ the central gas particle, be- 
cause the highest density region of the galaxy is likely to be 
the main contributor to both the UV luminosity and GRB/SN 
production. Note that a single gas particle in our simulation 
represents a gas mass of ^ 10^- IO^Mq, and /smooth.o is typi- 
cally on the order of kpc. 

In Figure[3] we show the distribution of simulated galaxies 
on the metallicity vs. extinction plane. Note that the metallic- 
ity used here is the mean gas metallicity of the whole galaxy 
(Zgaiaxy)- Only galaxics with My, > Kf -^M^ are shown in this 
figure, because the dust extinction in most of the lower mass 
galaxies are negligible {Eb-v ^ 0.01). The galaxies in the 
N400L34 run have lower Eb-v than those in the N400L100 
run for a given metallicity, indicating that ctism is lower in 
the N400 L34 run on average. Th is is consistent with the 
finding by Choi & Nagamine (2009) that the galaxies with a 
lower resolution run have higher gas mass fractions, because 
a higher resolution run can resolve higher density peaks and 
hence more gas is consumed for star formation. 

2.4. Luminosity Function of Simulated Galaxies 

Figure |4] shows the LFs of simulated galaxies in 
the N400L100 and N400L34 runs, comp uted using the 
GALA XEV population synthesis code ( Bruzual & Charloa 
120031) and the dust extinction model described in §|! 
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Each simulation can resolve dark matter halos and galaxies 
only in a limited mass range due to its limited box-size and 
resolution. In terms of the galaxy properties, the N400L100 
run contains a larger number of massive, brighter galaxies 
than the N400L34 run due to its larger box size. The N400L34 
run can resolve lower-mass, fainter galaxies better than the 
N400L100 run due to its higher resolution. 

The rest-frame Z?-band LPs are shown in the top left panel 
(dust extinction effect included) of Figure |4] with Poisson er- 
ror bars. It is clear that the N400L100 run misses the faint 
galaxies with Mg ^ -19, while the volume density of bright 
galaxies with Mg < -20 in the N400L34 run is smaller than 
that intheN400L100run. 

Since the two runs probe galaxies with different masses and 
luminosities, we combine the two LFs by taking the larger of 
the two runs at each Mb, and compare the result with the ob- 
servation by'Faber et al.' ("2007) in the top right panel of Fig- 
ure |4] The combined LF agrees quite well with the observa- 
tion when the dust extinction effect is included. 

We also compare the r est-frame UV LF with the observa- 
tion (iDahlen et al.l l2007h in the bottom panels of Figure H) 
Contrary to the case of rest-frame B-band LF, the discrepancy 
between the observation and simulation is significant at the 
bright-end of the rest-frame UV LF. We discuss the effects 
of this over pred iction on our luminosity PDF of GRB host 
galaxies in § 13.51 

Besides the bright-end overprediction, according to the 
comparison in FigureH] the simulation seems to overestimate 
the galaxy LF at any Muv by a factor of 1 .5 at z = 1 .0. How- 
ever, this over-estimation is probably within the range of un- 
certainties in our modeling and the cosmic variance in obser- 
vations. It should be noted that the vertical shift of LF (e.g., 
the uniform overprediction at any Muv) does not affect our 
luminosity PDF of GRB/SN host galaxies. 

2.5. Metallicity of Simulated Galaxies 

The distribution of simulated galaxies on the stellar mass 
vs. metallicity {Mi,-Z) plane is shown in the left panel of 
Figure|5] Here the metallicity means the mean gas metallicity 
of a galaxy. The N400L34 galaxies are distributed around 
the empirical formula at z = 1.0 proposed by Savagl io et ai.i 
12005), which agrees w e ll with the observational data points 
by Pe rez-Montero et all (120091) . The N400L100 galaxies are 
distributed on the more massive and luminous side than the 
N400L34 galaxies as expected from the luminosity functions. 
In both runs, more massive galaxies have higher metallicities 
than lower mass galaxies. 

We note that the N400L100 galaxies have somewhat lower 
metallicities on average than those in the N400L34 run for 
the same M^. This can be understood as follows. In our simu- 
lation, chemical enrichment is calculated as an instantaneous 
feedback from star formation, and the metallicity of a galaxy 
increases as the gas is turned into stars. Therefore the metal- 
licity of a galaxy is anti-correlated wi th th e gas mass fraction 
(/gas) of the galaxy. As mentioned in § 12.31 a higher resolution 
simulation can consume gas in star formation more efficiently 
than a lower resolution simulation. Furthermore, if we com- 
pare galaxies with similar M^,, galaxies with larger /gas have 
larger total (star and gas) mass, and thus reside in lager sim- 
ulation box. The N400L100 run has a higher f„^^ compared 
to the N400L34 run due to its larger box size and a lower 
resolution, leading to a lower metallicity of galaxies. As dis- 
cussed in § 12.41 the two simulation boxes complement each 
other, and each of the simulation cannot reproduce the obser- 



vations alone for a wide dynamic range. If we combine the 
two distributions of galaxies shown in the left panel of Fig- 
ure|5] the overall distribution of the simulated galaxies agrees 
with the observations. 

It should be noted that the mean gas metallicity (Zgaiaxy) of 
a simulated galaxy depends on the threshold density of gas 
particles to be included in galaxies, which is a parameter in 
the galaxy finding algorithm described in § 12.11 If we include 
lower density gas particles than the current density threshold, 
each simulated galaxy may have lower Zgaiaxy, because gas 
particles in the outskirts of a ga laxy have lower metallicity 
than those near the center (see § I2.2| |. However, the depen- 
dence of Zgaiaxy on the grouping density threshold does not 
affect our luminosity PDF of GRB/SN host galaxies. Low 
density gas particles are not forming stars, and hence they do 
not contribute to the GRB production and UV luminosity. 

The distribution of the simulated galaxies on the Luv - Z 
plane is also shown in the right panel of Figure |5] which is 
very similar to the M*-Z distribution. 

3. RESULTS 
3.1. Metallicity of GRB Sites 

In this section, we explore the metallicity of GRB sites. We 
first assume that the GRB occurrence rate in each simulated 
galaxy is proportional to the total SFR without a metallicity 
limit. Then we examine the models in which the GRB rate is 
proportional to the total SFR in gas particles with Z < Zcht. 
Our interest is in the luminosity PDF of GRB/SN host galax- 
ies, therefore we do not need to consider the absolute value of 
the event rate density. 

In order to investigate the effect of metallicity limit, we 
first examine the SFR-weighted distribution of gas particles 
as a function of metallicity. In Figure |6] we divide the galaxy 
population into three different samples according to the mean 
gas metallicity of the host galaxy (Zgaiaxy)- We can see that 
the metallicity of star forming gas particles has a wider range 
than that of Zgaiaxy, showing the large dispersion of gas metal- 
licity in each galaxy. We note that the high metallicity gas 
is concentrated in the high density regions of its host galaxy 
(see § 12. 2b where SFR is high. Therefore the peak metal- 
licity of SFR-weighted distribution is higher than the mean 
Zgaiaxy of the Sample. The distributions for both N400L34 and 
N400L100 runs are very similar to each other. 

Assuming that the GRB rate is proportional to the sum of 
SFR in gas particles with Z < Z^n, we can see in Figure |6] 
that, if we set Z^n = 0. 1 Z0, then the GRB event rate in galax- 
ies with Zgaiaxy ^ I.OZ0 (0.1 Zq) would be reduced to a few 
percent (40%) compared to the case of no Zdit. With Z^nt = 
0. 1 Zq, the event rate in galaxies with Zgaiaxy ^ 0.01 Zq would 
be almost unchanged from the case of no Zdit. The large fluc- 
tuation in the distribution for galaxies with Zgaiaxy ~ O.OIZq 
is due to the small number of gas particles included in such 
low-metallicity galaxies. 

3.2. Rest-frame UV Luminosity of the Host Galaxies 

Next, in order to discuss the luminosity PDF of host galax- 
ies, we consider the SFR-weighted fraction of galaxies as a 
function of rest-frame UV luminosity as shown in Figure |7] 
The luminosity PDF in the case of no Z^it is plotted in panel 
(a) for the N400L34 and N400L100 runs. As mentioned in 
§ 12.41 the two simulations have different resolving power of 
galaxies with different mass scales, therefore we combine the 
two distributions by taking the larger of the two runs in units 
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Fig. 4. — Luminosity functions of the simulated galaxies in the rest-frame B-band (top panels) and the rest-frame UV (2800 A, bottom panels). Top left panel: 
The solid (red) and dashed (blue) lines are for the N40()L34 and N400L100 runs, respectively. The error bars show the Poisson errors. Top right panel: The solid 
(magenta) line represents the combined LF of both the N400L34 and N400L100 runs (i.e., the larger of the two simulation at each magnitude). The LF without 
dust extinction (i.e., Eg-v = for all galaxies) is shown with the dotted (grey) line. The same line types are used i n the bottom panel s (rest-frame UV). The 
observed data points are taken from the DEEP2 survey for B-band iFaber et al. 2007), and the GOODS-S for the UV jPahlen et ani2007h . 



of Mq yr"' Mpc"^ mag"' as shown in panel (a). The distri- 
butions are normalized so that the integration of the combined 
luminosity PDF for each Z^^t is unity. Hereafter we discuss 
the combined luminosity PDFs unless otherwise mentioned. 

In panel (b) of Figure |7] the three combined luminosity 
PDFs of the simulated host galaxies with Zcm = none, Q.5Zq, 
and O.IZq are shown. The luminosity PDF shifts to the 
fainter magnitudes as Zcht becomes lower In the case of 
Zciit = O.IZq, the host galaxies are fainter than in the no-Zait 
case by ~ 1 magnitude. This result is consistent with the lumi- 
nosity PDFs of the observed GRB/SN host galaxies at z < 1.2 
(the F06 sample) as shown in panel (c). The redshift criterion 
of z < 1.2 was proposed in F06 to reduce the effect of galaxy 
evolution on the comparison between the GRB and SN host 
galaxies. The luminosity PDF of F06 was in the observed 
frame V-band, which is close to the UV in the galaxy's rest- 
frame. The median values of each distribution in panels (b) 
and (c) are indicated by the arrows. 

It should be noted that, using a larger observed sample, 
ISvensson et al.l (120101) found a smaller difference between the 
GRB and SN host galaxy luminosity PDFs than that claimed 
by F06. However, it is difficult to understand the selection 



effect in the GRB host sample of ISvensson et alj (120101) . be- 
cause a part of it is drawn from the GHostS project database, 
which consists of publicly available data of GRB host ob- 
servations done by various groups with differing instruments 
(ISavaglio et al.l 12009). On the other hand, the F06 sample 
is a complete sample of host galaxies for all 42 GRBs that 
have optical afterglow detections (> 3(j), therefore the sam- 
pling effect would be smaller in the F06 sample than in the 
Svensson et al. (2010) sample. Thus we use the GRB host 
galaxy sample by F06 as a reference data for the discus- 
sions in this paper. The 24 out of the 42 GRBs in the F06 
sample are at z < 1.2 . We note that the SN host galaxies 
in the Svensson et al.l (12010^ sample is significantly fainter 
than those in the F06 sample, and this is probably one of 
the causes for the smaller difference between the GRB and 
SN host galaxies. But the reaso n for the difference be tween 
the SN host samples of F06 and ISvensson et all ( 120101) is not 
known. 

Even the F06 sample would suffer from the observational 
bias against dusty host galaxies. We w ill disc uss the selection 
effect caused by the dust extinction in §14. 1 . II We note t hat the 
luminosity PDFs studied by F06 and ISvensson et al.l (120101) 
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Fig. 5. — Left: Distribution of galaxies on M^—Z plane. The galaxies in the N400L34 and N400L100 runs are shown with the solid (red) and dashed (blue) 
contours, respectively. The contours represent c/^;igai/(/y£/(logioM*)£/(logioZ) = 0.005, 0.05, and 0.5 [Ir' Mpc"^]. The stellar mass and metallicity of the observed 
galaxies at z ~ 1.0 fl'erez-Montero et al. 2009 ) and the empirical formula of Savaalio et al. 1 2005 ) for z = 1 .0 galaxies are plotted together The empirical formula 
of lSavaglio et alj i20b50 is based on the local and z = 0.7 galaxies with a mass range of ~ 10*^-10"Mq. Right: Same as the left panel, but on the Muy-Z 
plane, where M\]\ is the restframe UV magnitude. The contours represent iig^i/ dV dMi]\/d{logiQZ) = 0.002, 0.02, and 0.2 [h^ Mpc"^ mag"']. Colors are only 
for the online version. 




-3 -2-10 1 



Fig. 6. — SFR-weighted metallicity distributions of gas particles in the 
simulated galaxies with logio(Zgaiaxy/20) = 0.0 ±0.15, -1.0 it 0.15, and 
-2.0 ±0.15 are plotted with solid (red), dashed (green), and dotted (blue) 
lines, respectively. In our model, we assume that the GRB rate is propor- 
tional to the sum of SFR for the gas particles with Z < Z„n. The top and 
bottom panels are showing the distributions for the N400L34 and N400L100 
simulations. The large fluctuation in the particle metallicity distribution in 
galaxies with Zg^iaxy = 0.0 IZq (blue dotted line) is due to the small number 
of gas particles included in such low-metallicity galaxies. Colors are only for 
the online version. 

are at different wavelengths [rest-frame UV for F 06 versus 
rest-frame B & V-band for Svensson et al] (12010 )1. though 
this cannot explain the difference in their results. The dif- 
ference in the luminosity PDFs of GRB and SN host galax- 
ies should be larger in the optical than in UV, given that 
many of observed GR B host galaxies are extremely blue (e.g., 
iLe Floc'h et ani2003h . 

3.3. Cumulative Luminosity PDF 

We plot the cumulative luminosity PDFs of the simulated 
and observed host galaxies in Figure [8] The simulated PDFs 
are shown for the cases of Zcui/Zq = none, 0.5, 0.1, 0.025, 
& 0.001 (top to bottom). The observed host galaxy sample 



at z < 1.2 of F06 is plotted together as histograms. The shift 
of luminosity PDF to fainter magnitudes is clearly seen also 
in the cumulative plot, reproducing the difference in the ob- 
served GRB and SN-host galaxy luminosity PDFs. However, 
the distribution for no-Zdit case overpredicts the bright-end of 
the observed SN-host luminosity PDFs, and the Zait < O.IZ0 
case overpredicts the bright-end of the observed distribution 
of GRB host galaxies. This is expected from the overpre dic- 
tion of the UV LF at the bright-end as we discussed in § 12.41 
(see the bottom panel of Fig. |4]i. We discuss the effects of this 
overprediction on our conclusion in § 13.51 

3.4. Kolmogorov-Smimov Test 

To examine the consistency between the simulated and 
observed cumulative luminosity PDFs shown in Figure [8] 
we employ the one-sample Kolmogorov-Smirnov (K-S) test, 
which calculates the probability that a sample distribution is 
consistent with a given PDF. Here, the sample is the observed 
luminosity of GRB/SN host galaxies, and the PDFs are the 
results of our simulation. Figure |9] shows the K-S probability 
Pks (i e., one minus the rejection confidence level for consis- 
tency) as a function of Zait of simulated host galaxies. The 
Pks takes a value between and 1, and Pks < lO"'* means that 
the observed sample is dissimilar to the simulated one at a 3a 
level, while a larger Pks means a higher probability that the 
observed sample is derived from a PDF calculated in our sim- 
ulation. For the SN host galaxies, we find Pks > 0.05 when 
Zciit > O.O3Z0, and Pks > 0.2 when Z„it > 0.5Zq. On the 
other hand, for the GRB host galaxies, we find Pks > 0.05 
(0.2) when Zcdt/Z© < 0.1 (0.03). The KS probabiHties are 
the highest when Zcm/Zq = 0.005 and 0.5 for the GRB and 
SN host galaxies, respectively. This clearly supports the low- 
metallicity preference for GRB host galaxies. 

The luminosity PDF of the simulated host galaxies does not 
change significantly forZcrit < O.IZ0, and the agreement with 
the observed PDF is reasonably good for all Z^ni values for 
•Zcrit < O.IZq. Therefore we cannot put a lower limit to Z^ni 
based on this comparison. However our simulation cannot 
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Fig. 7. — (fl) Luminosity PDF of the simulated host galaxies (z = 1.0) with- 
out any metalHcity cutolf. The results of the N400L34 and N400L100 run are 
plotted with double-dotted (black) and double-dot dashed line (orange), re- 
spectively. The combined distribution of the two simulations (i.e., the larger 
of the two distributions) is plotted with soHd line (red), (b) Combined lu- 
minosity PDFs of the simulated host galaxies with no Zcit, Z^ii/Zq = 0.5, 
and 0.1 are plotted with the solid (red), dotted (magenta), and dot-dashed line 
(blue), respectively, (c) Luminosity PDFs of the observed sample of GRB 
host galaxies and SN host galaxies at z < 1.2 (F06) are plotted with dotted 
(dark blue) and dashed (dark red) histograms, respectively. The an'ows indi- 
cate the median of the distributions. Colors are only for online version. 
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Fig. 8. — Cumulative luminosity PDFs of the simulated host galaxies is 
shown together with those of the observed GRB/SN host galaxies of F06 at 
z < 1.2, based on the histograms shown in Fig. [7] The PDF without a metal- 
licity cutoff (no Zc,it) is shown with the solid line (red), and the those with 
Zcrit/Z0 = 0.5, 0.1, 0.025, and 0.001 are represented by the dot-dashed lines 
(blue) from top to bottom, respectively. Colors are only for online version. 
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Fig. 9. — K-S probabilities (i.e., 1.0 minus the rejection confidence level 
for consistency) of the luminosity PDFs of the observed and the simulated 
host galaxies is plotted against Z^rii. The symbols are for the SN host galaxies 
(red empty circles) and GRB host galaxies (blue filled circles). 

resolve the galaxies with M^, < IO^Mq, therefore it is possi- 
ble that the above result suffers from the resolution limit of 
our simulation. We should be cautious about the results with 
•Zcrit < O.IZ0, but our simulation contains galaxies as faint as 
Muv ~ -14 (see the bottom panels of Figure [Hi, which are 
well below the typical luminosity of the simulated host galax- 
ies with Zci-it < 0. 1Z0. 

3.5. Overprediction ofUV-bright Galaxies 

In §123] and Figure |4] we showed that our N400L100 sim- 
ulation overpredicts the bright-end of the rest-frame UV LF 
even after the dust extinction correction. We note however 
that our simulation agrees very well with the observed rest- 
frame B-band LF, which suggests that either our simulation 
overpredicts the amount of young stars compared to the older 
stellar population, or the two sets of observations are incon- 
sistent with each other, or else the Calzetti extinction curve 
that we adopted is not appropriate for the observed galaxies. 
The source of this discrepancy is unclear, therefore we focus 
only on the rest-frame UV LF in this subsection. 

To quantify the effect of the bright-end overprediction on 
our luminosity PDFs, we perform the same K-S test for the 
following two different test models that modify the shape of 
the rest-frame UV LF by hand. 

• Model 1: We double the gradient of the rest-frame UV 
LF at the bright-end (Muv < -20) as follows: 



Muv.model 1 =Muv-|-AMuv 



Muv < -20 
Muv > -20. 



(3) 



The SFR in galaxies with Muv < -20 is also modified 
in a consistent manner, motivated by the inference that the 
overprediction of the bright-end of the UV LF is caused by 
an overestimation of SFR in massive galaxies (e.g., due to 
lack of active galactic n ucleus feedback in our simulation; 
IChoi & Nagaminjl2010h . The change of SFR in each sim- 
ulated galaxy is as follows: 



A(logioSFR) = A(logioLuv) = - 



AMu' 
2.5 



(4) 



Model 2: We shift the entire UV LF uniformly by 
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Fig. 10. — Rest-frame UV luminosity functions of the simu lated galaxies 
in M odel 0, 1 & 2 are plotted together with the observation IDahlen et alj 
12007:) . The dust extinction model is included. The original result of our 
simulation (Model 0, identical to the combined simulation shown in the bot- 
tom right panel of Fig.|4) is shown with the dotted (black) line. The modi- 
fied LF with the doubled gradient of the bright-end slope (Model 1) and the 
horizontally shifted LF by AMuv = 10 (Model 2) are represented by solid 
(red) and dashed (blue) curves, respectively. Vertical dashed line represents 
Muv = ~20, which is the threshold of the bright-end modification in Model 
1. 

AMuv = 1 .0 to fit the observation. In this model, the luminos- 
ity PDF will also shift by the same amount of AMyv, without 
changing its shape. 

Both of these models are ad hoc modifications to fit the ob- 
served data points, but they would allow us to evaluate the ef- 
fects of overprediction of UV LF by our simulation. For com- 
parison purposes, we call the earlier calculations presented in 
§|33]-[l4]the Model 0. We show the LFs of Model 0, 1 & 
2 galaxies in Figure [TO] One sees that the overprediction of 
the bright-end in Model is removed in Model 1 & 2. The 
uniform overprediction at any Muv is not removed in Model 
1, but the uniform overprediction does not affect the luminos- 
it y PD F of GRB/SN host galaxies, as we have mentioned in 
§lMl 

Figure [TT] shows the K-S probabilities of Models 1 & 2 for 
the observed GRB/SN host galaxies. For the SN host galax- 
ies, the no-Zcrit case has the highest Pks, and Pks decreases 
with decreasing Zcrit. For the ORB host galaxies, the Model 1 
gives the highest probability of Pks ^ 0.6 at Zait = 0. 1 Zq, and 
Pks ~ 0.3-0.4 at Z^^ < O.IZq. In the case of Model 2, the 
trend is opposite of what we saw in Figure|9]for Zdit < 0.5 Z0 
in the sense that Pks increases with increasing Zcrit- But the 
Model 2 still prefers Zait = Q.5Zq to the no-Zait case. The 
Model 2 gives Pks < 0.2 across all Z^^i values for the SN host 
galaxies, which means that it doesn't agree well with the ob- 
servations of SN host galaxies compared to the Models and 
1. 

4. DISCUSSIONS 
4.1. Effect of Dust Extinction in the Host Galaxies 
4. 1 . 1 . Observational bias against dusty host galaxies 
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Fig. 1 1 . — Same as Fig.|9] but for the two different models of modified 
UV LF (Model 1 & 2). The Model 1 & 2 are represented by the squares and 
triangles, respectively. 




Fig. 12. — Same as Fig.[8] but here the galaxy is excluded from the sample 
if the average extinction value of galaxy exceeds £B-v,crit- We use Model 
in this figure. 

In this section we discuss the possibility that the systematic 
faintness and the low-metallicity of the observed GRB host 
galaxies is caused by the observational bias against dusty host 
galaxies, not by the intrinsic nature of GRB progenitors. It is 
difficult to identify the GRB host galaxy without a detection 
of an optical afterglow, and the afterglow flux may be reduced 
significantly by the dust extinction effect. In fact, some GRBs 
have optical to X-ray afterglow spectral index ^ox < 0.^ 
which is smaller than the expectation from t he standard ex- 
ternal shock model of GRB afterglows (e.g., Jakobsso n et alJ 
HoqI). We note that a part of the optically dai'k GRB could 
also be caused by the effects other than the dust extinction 
in their host galaxies, such as the attenuation b y neutral hy- 
drogen in the intergalactic rn edium at z > 6 (e.g.. lTotani et alj 
lIOOalNapamine et allllOOSl) . 

In order to examine the effect of dust extinction in the 
host galaxy (as opposed to the extinction in the immediate 
nearby environment of the GRB), here we exclude the galaxy 
from the sample of simulated host galaxies if Eb-v > EB-v.a-it, 

* The optical to X-ray afterglow spectral index /3ox is defined as fopt/fx = 
(Aopt/Ax)'^"*', where /opt and fx are the optical and X-ray fluxes, respec- 
tively. 
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Fig. 13 . — Same as Figs.l9land [m but here /^s is plotted against v,crit 
instead of Zcrii. 

wh ere E b-v is the average extinction of each galaxy computed 
in § 12.31 The GRB event rate in each galaxy is assumed to be 
proportional to the SFR in the galaxy. Figure[T2]shows the cu- 
mulative luminosity PDF of simulated host galaxies with var- 
ious Eb-yxxa- The no-ZiB-v.ci-it case (red solid line) represents 
the luminosity PDF of GRB host galaxies when the GRB rate 
simply traces the SFR, and it is identical to the no-Zcrit case in 
Figure [8] The luminosity PDF of the simulated host galaxies 
shifts to the fainter magnitudes with smaller Eb-v.cui, similarly 
to the Zcrit dependence shown in Figure [S] 

Figure [13] shows the K-S probability versus EB-v.mt for the 
luminosity PDF of the Model galaxies given in Figure [121 
together with the results of Model 1 & 2 galaxies. For Model 
& 1, a reasonable agreement (Pks > 0.2) between the sim- 
ulated and observed distributions of the SN host galaxies is 
found for a wide range of critical extinction {Eb-v.cui > 0.05), 
while for the GRB host galaxies the same level of agreement 
is achieved only when EB-v.mt < 0.03. 

The above requirement of Eb-vxA < 0.03 for the GRB host 
galaxies of the Models & 1 does not seem to agree with 
observations. Some follow up observations of optically dark 
GRBs suggest that a typical dust extinction to make a GRB 
optically dai'k is Ay - 1 (e.g.. lPerlev etani2009l Eb-v ^ 0.3), 
which is much larger than the above requirement of EB-v.ait < 
0.03. Such large values of critical extinction {EB-v.ciit > 0.1) 
do not significantly change the luminosity PDF of simulated 
host galaxies. A small dust extinction such as Eb-v < 0.03 
is unlikely to have a significant impact on the GRB optical 
afterglow observations. 

For the Model 2, the agreement with the observed GRB 
hosts is good (Pks = 0.48) at EB-v.ah = 0.1, however, f ks > 0.2 
is not obtained for the SN host galaxies with any values of 
EB-v.aii- The SFR-weighted luminosity PDF of the Model 2 
galaxies are too faint to agree with the observation of SN host 
galaxies even in the no-EB-v.ait case. Hence the difference of 
luminosity PDFs between the GRB host and SN host galaxies 
cannot be reproduced by the Model 2. 

A caveat here is that the Eb-v is computed in our model 
assuming simplified structure of ISM in each galaxy without 
considering smaller scale structures such as molecular clouds. 
Though our model of dust extinction reproduces the observed 
B-band LF of field galaxies, the model may be not representa- 
tive of dust extinction for lines of sight to GRBs, if small scale 
structures near GRBs play an important role in the extinction 
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Fig. 14. — Fraction of optically dai'k GRBs at z= 1.0 in our simulation (i.e., 
the fraction of GRB events in galaxies with Eb-v > ^B-v.crit) as a function 
of Sfi-vciit- The results for three different values of Z(.,it are shown. The 
observed range of fractions for the GRBs at all redshifts (Cenko et al. 200^; 
IZheng et alJ2009l) is indicated with the shade. We use Model in this figure. 

of GRB optical afterglows. Recent observations of optically 
dark GRBs actually sugg est that the effect of dust on small 
scales may be importan t (jRol et al.ll2007| ICastro-Tirado et aD 
2007; Perlev et alj2009H Hashimoto et al. 2010; Holland et al] 
2010). If the typical extinction in lines of sight to GRBs is 
greater than the average extinction in its host galaxy by one 
order of magnitude but yet proportional to the average extinc- 
tion, the luminosity PDF of GRB host galaxies possibly be 
reproduced by the observational bias against dusty host galax- 
ies. 

In this section, we have considered the possibility that the 
observational bias against dusty host galaxies is causing the 
systematic faintness of GRB host galaxies. However, we note 
that the SN observations possibly suffer more from dust ex- 
tinction than the optical GRB afterglow observations. In that 
case, the observational bias against dusty host galaxies would 
make the SN host galaxy sample fainter than the GRB host 
galaxies, which is opposite to the observational result. Our 
simulation does not contradict with the observations even if 
the SN observations suffer from dust extinction, because the 
simulated luminosity PDF with EB-v,cnt > 0.1 agrees with 
the SN host observations reasonably well (Pks > 0.2) for the 
Model 0. 

4. 1.2. Fraction of Optically Dark GRBs 

The physical nature of optica lly dark GRBs is al so an im- 
portant issue. A recent study by Cenko et al. ('2009) reported 
that the fraction of optically dark GRBs in the entire observed 
GRB sample is - 50%, while Zheng et alJ (12009) found that 
the fraction is ^ 10-20%. Note that these fractions are based 
on the entire GRB sample that span a wide redshift range, be- 
cause it is difficult to constrain the redshift of optically dark 
GRBs. 

In Figure [141 we show the fraction of GRB event rate in 
galaxies with Eb-v > Eb-vxA at z = 1.0 in our simulation 
for three different values of Za-it- Assuming that a GRB 
at z = 1.0 becomes optically d ark when the host galaxy has 
Eb-v > EB-v.crit = 0.3 (Ay > l. lPerlev et al1l2009l) . our simu- 
lation predicts that the fraction of optically dark GRBs is < 
10% for all cases of Zan- For the no-Zoit case (red solid Une), 
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the dark fraction is close to 10% for EB-v,cnt = .3, and is al- 
most c onsistent with the observational result of IZheng et alj 
(120091) . However, if Z„it < 0.5Zq (as suggested in §[3]l, our 
simulation predicts that the fraction of optically dark GRBs 
is < 5% for Eg-v.cnt = 0.3, which is much smaller than the 
observed fraction. 

The small fraction of optically dark GRBs in our simula- 
tion compared to the observations can be explained in two 
different ways. The first possibility is that the optically dark 
GRBs mainly suffer from the absorption in the optical wave- 
length due to dust associated with small sc ale st ructures, not 
from the galaxy scale dust as modeled in § 12.31 The second 
possibility is that many of the observed optically dark GRBs 
originate at different redshifts other than z 1 .0, although the 
redshift distribution of GRBs with known redshifts that are 
typically optically non-dark, peaks at z ^ 1 .0. The fractional 
contribution of dusty galaxies to the cosmic SFR would be 
different at different redshifts. Furthermore, the afterglow of 
a GRB at z > 6 would be attenuated by neutral intergalactic 
medium in the observer-frame optical bands. 

It is unclear which of the above two scenarios might 
be correct at this point. Som e observational studies 
JCfistro-Tirado et al. 20071 iPerley et al. 2009; Holland et al. 
120 lOi) find that the host galaxies of optically dark GRBs have 
bluer colors than those expected from the dust reddening es- 
timated from the afterglow spectra. On the other hand, the 
host galaxies that is globally dusty are also found for several 
optica lly dark GRBs (Levesque et al. 2010; Has himoto et alj 
l2O70h . 

4.2. Comparison with Previous Works 

Several studies have been performed to interpret the ob- 
served luminosity PDFs of GRB and SN host galaxies. In 
this section we compare our study with previous results. 

Wolf & Podsiadlowski ( 2007) is one of the important stud- 
ies in this context. They reproduced the luminosity PDFs 
of the observed GRB and SN host galaxies based on the 
empirical relations of galaxy properties at z ~ 0.7, such 
as the luminosity-metallicity relationship. In their study, 
the model in which the GRBs occur in high metallicity re- 
gions (Z ^ Zq) reproduces the luminosity PDF of GRB 
host galaxy the best, and the models which requires GRBs 
to occur only in sub-solar metallicity regions do not repro- 
duce the observations. They assumed a power-law relation 
for the luminosity-metallicity relation of galaxies following 
iKobulnicky & Kewleyl (120041) . However this relationship is 
not well constrained by current observations at redshifts of 
our concern, and in the work of Kobulnicky & Kewley ( 2004), 
different fitting methods gave different set of best fit param- 
eters for a set of observed data, suggesting large uncertain- 
ties in the fit. We also note that the results probably depend 
on the assumed dispersion around the power-law luminosity- 
metallicity relatio n. 

Kocevski et alJ (l2009h studied the stellar mass PDF 

of GRB host galax ies in a similar way to that of 

IWolf & Podsiadlowskil (l2007h and argued for Z^nt > 0.5 Zg. 
The mass-metalUcity relationship of galaxies a re better under- 
stood than the luminosity-metallicity relation (iTremonti et alJ 
|2()04; Savaglio et al. 2005). However, it is difficult to un- 
derstand the selection effects of the currently available sam- 
ples o f GRB host galaxies with known M*. iKocevski et alj 
(120091) discussed the observed sa mples of GRB host galax- 
ies with known M^, collected by Castro Ceron et al] ( l2008h 
and ISavagUo et al, (.2009.) . Both of the samples were col- 



lected from publicl y availa ble data of GRB host observations. 
ICastro Ceron et alj ( 120081) collected all GRB host galaxies 
with rest-fram e /T-band data including upper limits, while 
ISavaglio et al.l (|2b09) selected their sample requiring multi- 
band detection of the host galaxies. The observations of GRB 
host galaxies in the samples are originally done by various 
groups with differing instrument, and hence the selection ef- 
fects in the samples are hardly understood. We also note that 
the constraints on the stellar masses for some of the GRB host 
galaxies in their sample are not accurate enough to quantita- 
tively discuss the low-metall icity preference. 

The approaches o f IWoTf & Podsiadlowskil (|2007!) and 
IKocevski et al.l (120091) are entirely based on the empirical re- 
lations, and are fundamentally different from our approach 
based on numerical simulations. As discussed above, one 
difficulty in studying the luminosity PDFs of GRB/SN host 
galaxies is that properties of galaxies are not observation- 
ally well understood at redshifts where the observed sample 
is available. Numerical simulation approach is powerful in 
this context, because it gives properties of galaxies such as 
luminosity, metallicity, and SFR without assuming uncertain 
relations between them. 

Some numerical studies have inv estiga ted the lumi nosity 
of GRB host galax ies (iNuza et alJ 120071: iLapi et alJ 120081: 
ICampisi et al.l 120091) . However, quantitative comparisons of 
the predicted luminosity PDF of GRB host galaxies and the 
observe d luminosity PDFs have not been performed in these 
studies. iNuza et aP (f2007) did not reproduce the luminosity 
PDFs of both GRB and SN hosts, presumably o wing to the 
small box size of their s i mulati on (lO/z"' Mpc). iLapi et al] 
(2008) and ICampisi et all (l2009h have shown that the fainter 
galaxies have systematically lower met allicities in the i r mod - 
els. However, in Lapi et al. {200^ and lCampisi et alj (120091) . 
the luminosity PDFs are calculated simply by selecting galax- 
ies that contain young low-metallicity stars without taking the 
difference of GRB rate among galaxies into account, therefore 
the predicted luminosity PDFs cannot be directly compared 
with the observations of GRB host galaxies. 

In this context, our present work is unique in the sense 
that it can quantitatively reproduce the observations of both 
GRB and SN host galaxies using a self-consistent cosmolog- 
ical hydrodynamic simulations of galaxy formation. This has 
been enabled by applying a widely used GRB rate model to 
our cosmological simulations with large box sizes, which can 
properly compute the properties of galaxies at z 1 . 

5. CONCLUSIONS 

Using cosmological SPH simulations, we have examined 
the relation between the metallicity dependence of GRBs 
and the luminosity PDF of GRB host galaxies. Our results 
suggest that the observed difference in the UV luminosity 
PDFs of GRB/SN host galaxies can be explained by the low- 
metallicity preference of GRBs. 

We find that in our simulation, the luminosity PDF of GRB 
host galaxies agree with the observed one when Z^nt < Q.SZq, 
while the observed luminosity PDF of SN host galaxies can 
be reproduced without a metallicity dependence of SN rate. 
The suggested value of Z„\t for the GRBs by our study is 
consistent with the suggestion from stellar evolution models 
(lYoon & Langer 2005; Woosley & Heger 2006), contrary to 
the results of previo us studies ( Wolf & Podsiadlowskil 120071 : 
IKocevski et alJ2009l) . Though our simulation suffers from the 
seeming overprediction of the UV luminosity function at z = 1 
at the bright end, we have explicitly demonstrated that the 
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quantitative agreement between the simulated and observed 
luminosity PDFs is not a product of this seeming overesti- 
mate of the UV-bright galaxies in our simulation. However, 
we should be cautious about our results when Zcrit < O.IZq, 
since they may be affected by the resolution limit of our sim- 
ulation. 

We have also discussed the effect of observational bias 
against the host galaxies of optically dark GRB on the lu- 
minosity PDFs. The dust bias causes a shift of luminosity 
PDFs toward fainter magnitudes, similarly to the effect of the 
low-metallicity preference of GRBs. However, to obtain an 
acceptable fit to the observed luminosity PDF of GRB host 
galaxies, we had to assume that an optical afterglow becomes 
significantly faint when Eg-v ^ 0.03 in the host galaxy. Such 
a small extinction is unlikely to cause a significant effect for 
the optical follow-up observations of GRBs, and therefore it 
would be difficult to explain the observed luminosity PDFs of 
GRB and SN host galaxies only by the overall dust bias of the 
host galaxies. Nevertheless it is possible that our simplified 
model of dust extinction is not sufficient to reproduce the dust 
extinction of GRB optical afterglows. 

Assuming that a GRB becomes optically dark when the ex- 
tinction of its h ost galaxy is Eb-v ^ 0.3 as suggested by the 
observation of IPerley et aH ([2009), our simulation predicts 
that the fraction of optically dark GRBs at z = 1.0 is < 5% 
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